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I. INTRODUCTION 

The recent observations of neutrino oscillations [JQ-El 
and reported, although not confirmed, discovery of neu- 
trinoless double beta decay Q give strong motivation for 
studying some aspects of physics beyond the Standard 
Model (SM). One of the most popular approach to such 
a physics is the incorporation of supersymmetry, which 
leads, in the simplest case, to the Minimal Supersym- 
metric Standard Model (MSSM) 0. This framework is 
very attractive for theorists since it is the simplest su- 
persymmetric model possible that leads to unification of 
the coupling constants. It also provides a natural way to 
incorporate of lepton and baryon number violating pro- 
cesses, which is impossible in SM. 

It was widely approved to neglect the terms violating 
lepton and baryon numbers in actual calculations, due 
to smallness of the relevant coupling constants 0, Q, @ • 
However, even if perfectly valid in low energies, this as- 
sumption could have lead to errorneus predictions in high 
energy limits. We show by explicit simulations that this 
approach is justified in most cases. 

The paper is organized as follows. In the next section 
we present the model. We also explain our conventions 
and notations here. In Section III the results are pre- 
sented and discussed. A short conclusion is given at the 
end. 



II. THEORY 

The minimal (in the particle and interaction content) 
supersymmetric extension of SM means introducing a su- 
perpartner for each particle and adding a second Higgs 
doublet (together with its higgsinos). 

We closely follow the conventions used in [jj . The 
MSSM superpotential can be written in the form 



W 



MSSM 



e a b[(Y E)ijI j iH 1 Ej + (Y D )ijQf x HiDj X 
(Y u ) ij Qt*HlU jx +iiHiH»\, (1) 
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where Y's are 3x3 Yukawa matrices. L and Q are the 
SU(2) left-handed doublets while E, U and D denote the 
right-handed lepton, up-quark and down-quark SU(2) 
singlets, respectively. Hi and H2 mean two Higgs dou- 
blets. We have introduced color indices x,y,z — 1,2,3, 
generation indices i,j,k — 1,2,3 and the SU(2) spinor 
indices a,b,c = 1, 2. 

Since supersymmetry is not observed, we need to intro- 
duce some mechanism to break it. The simplest possibil- 
ity is the so-called soft breaking given by the Lagrangian 



[•soft 



£ab[{J±E)ij 



lth\i 



{A D ) ljq rh\d 3X 



(Au^q^hiu^ + Bfihthi], 



(2) 



where lowercase letters stand for scalar components of 
respective chiral superfields, and 3x3 matrices A as well 
as Bfi are the soft breaking coupling constants. We in- 
troduce also a scalar mass term of the form 



• 2 .h\h 2 + q'mgq 



Hi 



-,t r 



•|t m £z 



ij/" TuuijjU' T eiiige . (3) 

Looking at (JTJ it is clearly visible, that W MSSM pre- 
serves a discrete symmetry called i?-parity and defined 
as 

Rp = (_l)3B+^+2S ; 

where B, L, and S are the baryon, lepton and spin num- 
bers, respectively. From the definition follows that or- 
dinary particles have R p = +1 whereas supersymmetric 
(SUSY) particles have Rp — — 1 . The consequence of Rp 
conservation is a stable lightest SUSY particle, which can 
be a serious candidate for cold dark matter, and a stable 
proton. 

However, nothing motivates theoretically R p conserva- 
tion, especially that in order to avoid rapid proton de- 
cay, at least the lepton number parity L p = ( — i) L + 2S 
or the baryon number parity B p — (— \^ 3B + 2S must be 
conserved, but not necessarily both. It follows that we 
can add new terms to the superpotential, which violate 
separately the lepton and baryon number: 
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tab 2^ Ek ^ i ^i^ l> 3^ k (A-D k )ijLiQ* b Dki 

^(Au^UfDpt + e^LfHl (4) 
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The full superpotential becomes now 

W = W MSSM + W n " (5) 

and defines a new model, often called $ p -MSSM. It in- 
troduces nine new Yukawa coupling constants but, as 
mentioned above, proton stability requires to reject ei- 
ther lepton number or baryon number violating terms. If 
one wants to apply this model to the description of the 
neutrinoless double beta decay 

A(Z, N) -> A(Z + 2, N - 2) + 2e' , 

in which the lepton number is violated by two units, the 
matrices Ajji have to be set equal to zero. In this paper 
we will not do that in order to stay completely general 
and to be able to count contributions coming from all 
sources. 

All the parameters (coupling constants and masses) of 
^?p-MSSM are running parameters in the sense of renor- 
malization group, that is they obey a set of equations of 
the form 

(n) 

where t — log E and E is the renormalization scale; /3i 
is the n-loop renormalization group beta function for the 
parameter x. Our previous calculations |g, llll Il2j were 
based on 1-loop MSSM renormalization group equations 
because we roughly estimated that the corrections com- 
ing from i?-parity violation are of the order of few per- 
cent. We check this assumption in detail in the next 
section. 



III. RESULTS 

The relevant 1- and 2-loop equations for MSSM can be 
found in with$ p -MSSM corrections in [lO]. Our pro- 
cedure of runnig the parameters consists of few steps. We 
take into account mass thresholds where SUSY particles 
start to contribute [lj| and use the SM 2-loop renormal- 
ization group equations (RGE) below appropriate thresh- 
old and 2-loop |? p -MSSM RGE above it. Initially all the 
thresholds are set to 1 TeV and are dynamically modified 
during the running of mass parameters. After evolving 
the dimensionless couplings from the electroweak scale 
tyl z up to mcuT ~ 10 16 GeV, we unify the masses of 
gauginos, sfermions and squarks to be m Q — 500 GeV. 
The trilinear soft couplings at ttigut are set according 
to formula 

Ai = AqYi, 

with Ao — 500. Next, all the quantities are evolved down 
to mz- The numbers used in simulation serve as an ex- 
ample. Our investigation showed, that the same conclu- 
sions can be drawn for a wide range of mn and Aq. 



TABLE I: Mass matrices of supersymmetric particles (third 
family) and values of the coupling constants, calculated using 
1-loop MSSM and 2-loop $ P -MSSM renormalization group 
equations. All supersymmetric masses are unified (500 GeV) 
at the GUT scale 10 16 GeV. Other MSSM parameters were 
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The values of Yukawa couplings Y at mz are given by 
lepton and quark mass matrices M [§J 

M v - <#°)S^Y£S^, 

M D = (HI)S Dr YI& Dl , (7) 

M E = {H° x )S Er Y t e & El , 

where (H®) and (H®) are the neutral Higgs vacuum ex- 
pectation values. S matrices perform diagonalization so 
that one obtains eigenstates in the mass representation. 
However, we had to choose some low energy values for 
the remaining i?-parity violating couplings A. There is 
no hint concerning their actual numerical values except 
that they should be small, since lepton number violating 
processes are very rare. We checked several possibilities 
and found, that in order to obtain physically reliable re- 
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FIG. 1: _R-parity violating Yukawa couplings evolution. All were set to the same value at the mz scale and evaluated using 
the 2-loop ^p-MSSM RGE to the m GU T scale. 



suits, all these couplings should obey 

(A) u < 0.015, (*= 1,2,3), 



(8) 



otherwise high vacuum expectation values are generated 
from the down quark mass evolution. What is more, val- 
ues of other parameters nearly did not change for various 
A's. From we can make use of the following relation 



(A D i) n <£ 



10 2 



-1/4 



y 



t, 



Qv2(3 — exp 
1/2 



(9) 



Using the £ parameter for 76 Ge isotope for the unification 
mass 500 GeV, and inserting the experimentally known 
half-life given by the Heidelberg-Moscow collaboration 
1 T ^ iHMexp) = (0.8 - 18.3) x 10 25 2/, we obtain 



(A m ) u < 0.275, 



(10) 



which is much weaker than our constraintmit (8). If one 
uses limit (8) the half-life to be of the order of ~ 1.3 x 
10 30 y which is rahter an overestimation. On the other 
hand, the newest data from the IGEX Collaboration 0] 
sets the lower bound on the half-life in germanium to 

be T OjmiGEXe X p) > x 1Q25y withQut bQunds from 

above. 



Values for various parameters in both 1-loop MSSM 
and 2-loop $ p -MSSM cases are given in Tab|l| The 
SU(3) C x SU(2)l x U{1)y dimensionless coupling con- 
stants gi are represented here in the usual way by ai — 
gf/Air. Whenever the parameter has a matrix form, it 
has been diagonalized and the entry {33} is listed in 
the table. Masses are given respectively for gauginos, 
left-handed Higgs doublets, left-handed squark doublets, 
right-handed up- and down-type squark singlets, left- 
handed slepton doublets, and right-handed slepton sin- 
glets, with respect to the SU(2) group. 

In general, as expected, all parameters are changed 
by at most few percent. One should keep in mind that 
the new results include two corrections: the 2-loop terms 
and the terms connected with i?-parity violation. Our in- 
vestigation shows, that the differences are mainly driven 
by the 2-loop corrections, whereas the i?-parity terms 
change the results on the third decimal place. 

This situation can be understood from the evolution 
of -Rp-violating couplings. In Fig^ an example of run- 
ning of A coupling constants is shown. It is essential, 
that their values decrease with increasing energy. As a 
consequence the contributions coming from A's can be 
safely neglected in actual calculations. It is, however, a 
surprising result, since normally one would expect to see 
exotic processes at high energies. Fig^ clearly shows, 
that to observe a lepton number violating process, like 
the neutrinoless double beta decay, one should search for 
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it in low temperatures rather than in accelerators. 
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In conclusion, we have presented the spectrum of su- 
persymmetric particles and coupling constants using the 
1-loop MSSM and 2-loop ^ p -MSSM renormalization 
group equations. Inclusion of the i?-parity violating 
terms does not change the results significantly. On the 
other hand we confirmed the well known rule that the 2- 
loop corrections should be used in all cases, where quan- 
titative rather than qualitative features are important. 
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